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1. INTRODUCTION TO COMPLEX QUANTUM ORDERS 

One of the most challenging and interesting phenomena in modern condensed 
matter physics is the one emerging from strongly coupled systems as a result of com- 
peting interactions. Indeed, the multiplicity of distinct and novel quantum phases 
observed experimentally confront us with new paradigms that question our under- 
standing of the fundamental principles behind such complex phenomena. For exam- 
ple, whether the mechanism controlling the coexistence and/or competition between 
magnetism and superconductivity or Bose-Einstein condensation has the same phys- 
ical origin in different classes of materials is still an open question [1]. It is believed 
that the physics in these materials is strongly influenced by their proximity to quan- 
tum phase transitions and, in particular, to quantum criticality. A quantum phase 
transition is characterized by the qualitative changes of the macroscopic state of the 
system induced by tuning parameters of its Hamiltonian. On general grounds, the 
very notion and nature of entanglement is at the core of the problem [2]. 

From the theoretical viewpoint the hurdle is in the presence of non-linear cou- 
plings, non-perturbative phenomena, and a panoply of competing quantum orders. 
These systems happened to be strongly correlated since no obvious small coupling 
constant exists, and consequently exhibit high sensitivity to small parameter changes. 
It is then clear the importance of developing a methodology, based on qualitatively 
new concepts, that treats all possible competing orders on an equal footing with no 
privileged fixed-point phenomenon. Despite great advances there is a lack of a sys- 
tematic and reliable methodology to study and predict the behavior of these complex 
systems. It is a purpose of this work to present a promising step in that direction. 

In the quantum description of matter, a physical system is naturally associated 
with a language of operators. We have previously developed an algebraic framework 
for interacting quantum systems that let us study complex phenomena characterized 
by the coexistence and competition of various broken symmetry states [3-5], and 
proved a theorem that allowed us to connect all possible languages used in the quan- 
tum description of matter. Connecting the various languages through isomorphic 



mappings enable us to relate seemingly different physical phenomena, unveil hidden 
symmetries (i.e., uncover the accidental degeneracies of the original physical system), 
and, in some limiting cases, obtain the exact spectrum of the problem (or of a set 
of orthogonal subspaces). The ultimate goal was to use that framework to explore 
those unconventional complex states of matter from a unified perspective. Given the 
space limitations we cannot review these concepts but refer the reader to a recent 
review article on the subject [5]. 

The modern theory of phase transitions starts with Landau's pioneering work 
in 1937 [6]. One of his achievements was the realization of the fundamental relation 
between spontaneous symmetry breaking and the order parameter (OP) that mea- 
sures this violation, thus giving simple prescriptions to describe order in terms of 
irreducible representations of the symmetry group involved. Another was the devel- 
opment of a phenomenological calculational scheme to study the behavior of systems 
near a phase transition. Landau's theory has been successfully applied to study phase 
transitions where thermal fluctuations are most relevant. Certainly, the theory was 
not designed to study zero-temperature (quantum) phase transitions. 

In previous work [3,4] we outlined a framework to identify OPs based upon 
isomorphic mappings to a hierarchical language (HL) defined by the set of operators 
which in the fundamental representation (of dimension D) has the largest number of 
symmetry generators of the group. Any local operator can be expressed as a linear 
combination of the generators of the HL. The building of the HL depends upon the 
dimension D of the local Hilbert space, Tij, modeling the physical phenomena. For 
instance, if one is modeling a doped antiferromagnetic (AF) insulator with a t-J 
Hamiltonian [7], then D=3 (i.e., there are three possible states per site) and a HL 
is generated by a basis of the Lie algebra su{3) in the fundamental representation 
[3,4]. As explained and proved in Refs. [3-5], there is always a HL associated to each 
physical problem. These ideas complement Landau's concept of an OP providing a 
mechanism to reveal them, something that is outside the groundwork of his theory. 
Indeed, Landau's theory does not say what the OPs should be in a general situation. 

As mentioned above, these isomorphic mappings not only unveil hidden sym- 
metries of the original physical system but also manifestly establish equivalences 
between seemingly unrelated physical phenomena. Nonetheless, this is not sufficient 
to determine the exact phase diagram of the problem: One has to resort to either 
numerical simulations with their well-known limitations or, as will be shown in the 
present paper, to a guided approximation which at least preserves the qualitative 
nature of the possible thermodynamic states. A key observation in this regard is the 
fact that typical model Hamiltonian operators written in the HL become quadratic 
in the symmetry generators of the hierarchical group, and this result is independent 
of the group of symmetries of the Hamiltonian. 

This latter result suggests a simple approximation, based upon group theoretical 
grounds, which deals with competing orders on an equal footing and will be termed 
hierarchical mean-field theory (HMFT). The HMFT is distinctly suitable when the 
various phases displayed by a system are the result of competing interactions and 
non-linear couplings of their constituents matter fields. 



2. HIERARCHICAL MEAN-FIELD THEORY AT WORK 



In the rest of the paper we will expand on a technique to build approximate phase 
diagrams, dubbed HMFT. To avoid excess of tedious formalism we will describe 
the methodology by example. We have chosen two representative examples. The 
first displays two quantum phases separated by a bosonic metal-insulator transition 
induced by particle interactions. The second shows an AF, a superfluid and, in 
addition, a coexistence phase as a function of the particle density and temperature. 

2.1 A Superfluid-Insulator Transition 



Let us determine the zero temperature phase diagram of a very simple model 
displaying a Mott insulating to superfluid transition. The model we refer to is a 
modification of the Bose-Hubbard Hamiltonian [8] {U, t > 0) 

H = -'^tY.^aUi + gl9i) + U ^(nj - n + l)(nj - n) + ^ ej nj , 

(i,j> j j 

= -2^E(^?^j + alOi) + t^E"j("J - 1) + E/^J + NsUh{n - 1) , (1) 

(i,j> j j 

where (i, j) stands for nearest-neighbor sites (bond) in an otherwise regular A^g-sites 
lattice of coordination z, //j = 2U (1 — n) -t-ej, and the ^-particles, instead of canonical 
bosons, represent bosons satisfying the following commutation relations ((^j)^ = 

[9-v fj ] = '^ij ~ "-j) ' [^i' = <^ij 9} , 

defining an algebra, with n > 1 a positive integer and the nilpotency condition 
(S'j )^^^ = (gj)^ — meaning that one can accommodate up to p = 2 particles 
per mode j. This indicates that the local Hilbert space Tij is three-dimensional. A 
possible basis is {|n — 1), |n), \n -\- 1)}, and the operators acting on this basis 



fifjjn — 1) = , gj\n — 1) = \/n |n) , n^\n — 1) — {n — l)\n — 1), 

g^\n) = \ffi |n — 1) , g\\n) = \ffi |n -|- 1) , = n\n), (3) 

5rj|n+ 1) = |n) , ^j|n + l) = 0, nj|n + 1) = (n + l)|n + 1), 

behave as harmonic-oscillator-like operators acting on a three-dimensional space. 
Since D = 3 one can easily determine the isomorphic mapping between the (7-particle 
and the su(2) algebra in the S = \ representation. The transformation is given by 

' 9\ = , 

I o. = Sr (5) 
nj-n = S? , 

while the Hamiltonian operator in the spin language reads 

H = -tnJ2 {s^Sr + SrS+) + U ^(5?)^ + ^(C/ + ej) S? . (6) 

(i,j> j j 



Before translating the problem into the most fundamental language, i.e., the HL, 
let us make a detour and re-express H using the spin-1/2 (cr =t, [) Jordan- Wigner 
bosons [9], 

St ^ {h\^ +\{) , 



H = -2th 



where the number operator nj 
by these hard-core bosons is [5]: ^\^-,'j-^^t 



j j 

^jT + ^jJ. (^jo- = ^jcr^jcr)' algebra satisfied 

ha'] = 0, K,bU = - 2^j- - ^j^) (if 

As explained in the introduction, the first step in determining its phase diagram 
consists of re-writing in a HL. The latter is realized by S't/(3)-spin generators in 
the fundamental representation, and its mapping to the hard-core boson language 
can be compactly written as [3] 
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The three components s? — (nji- — njj^)/2, Sj" 
the spin su{2) subalgebra, i.e., they are the components of the local magnetization. 
The five additional components correspond to the Bose-Einstein condensate and the 
charge density wave local OPs. In the HL, H represents a Heisenberg-like Hamilto- 
nian [10] in the presence of an external magnetic field (J^^, = 7,^^) 

(i,j> j j 

(10) 

with Jqq = Jii = J22 = J12 = 0, Joi = J02 = —2tn. Note that through the mapping 
we transformed an interacting problem into another problem that is quadratic in 
the basis of the algebra su{3) in the fundamental representation, but which is not 
necessarily SU{3) symmetric. The idea behind the HMFT is to perform an approxi- 
mation which deals with all possible local OPs on an equal footing with no privileged 
symmetry axes and, hopefully, retains the qualitative topology of the phase diagram. 
With H written in the HL one immediately realizes that the simplest HMFT can 
be achieved if we re- write H in terms of SU (3) Schwinger- Wigner (SW) bosons (3 
flavors CK =1, 0, t) [7]. The mapping is expressed as 
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with the SW bosons bj^ satisfying the constraint riji + njo + nj-|- = 1. In this way 
iy=-2tnE (:4yA^y:-Si.jS,ij)-C7j]njo + J](C/ + ej)(njT-nji) , (12) 

{i,i),cr j j 

where : : denotes normal ordering and 



4ii = blbl - blobl 



(13) 



If ej is translationally invariant, then = S? = X]j(nj-|- — n^i) is a constant of 
motion {[H, S^] =0). In the following we will only consider the case 5'^ = 0. 

Since the su{N) languages provide a complete set of HLs [11], any model Hamil- 
tonian can be written in a similar fashion once we identify the appropriate HL and 
apply the corresponding SW mapping in the fundamental representation (the order- 
ing operators will, of course, have a different meaning and algebraic expressions). 
The key point is that the Hamiltonian operator in the HL becomes quadratic in the 
symmetry generators of the hierarchical group {SU{3) in the present case). 

The idea behind any MF approximation is to disentangle interaction terms into 
quadratic ones replacing some of the elementary mode operators by their mean value. 
The crux of our HMFT is that the approximation is done in the HL where all possible 
local OPs are treated on an equal footing and the number of operators replaced by 
their mean value is minimized since the Hamiltonian is quadratic in the symmetry 
generators. In this way, the information required is minimal. In mathematical terms, 
given Oj. o„ ^ (ot.^o„ + ol.{0,^) - + (Oj, - {Ol.)){0,^ - (O^)), for an 

arbitrary bond-operator , the approximation amounts to neglect the latter fluctu- 
ations, i.e., ClCy ^ (Ojj)Oij + - (C»it)(Cij). An important result is that aU 
local OPs are equally treated and, moreover, symmetries of the original Hamiltonian 
related to the OPs are not broken explicitly in certain limits. 

The resulting MF Hamiltonian together with the SW-boson constraint (with 
Lagrange multiplier X) H = Hmf + -^X]j(^jJ. + ^jo + ^jt) ^^ads [10] 

H = -2tn [A{Al. + A^,^ - iB{Bl,. - B^,^)] - C/ ^ njo + A ^ nj^ 

= Y. ^ks^ko + As &L^Lko + H.C.] [(A - U)n^o +A(nks + nka)],(14) 
keBZ keBZ 

where the sum of momenta k is performed over the first Brillouin zone (BZ), ^>j^g(a) ~ 

;^[6|,| ± 6]^J, nica = ^ka^ka (^ka'^ represent Fourier transformed modes). = 

-2^/2 tnA7k, Kb = 2^2 tnB%, with 7k = 2 Xl)!=i cos(A;^) and 7k = 2 X))l=i sin(A;^). 
The resulting self-consistent MF equations are 

\/2 
* keBZ 

^ = E 7k(&L^Lko)MF , (15) 
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1 = ^ E E<'^ka)MF 
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On the other hand, it can be shown that the case with vanishing S-ordering is a 
stable solution. In the zero-temperature limit all particles condense into the k = 
mode. This condensation corresponds to the appearance of a superfluid phase of 
5f-particles (or XY S=l ferromagnetism in the spin language) and it is manifested 
in the non-zero value of A: 

This expression can be easily obtained from Eq. (15) considering that at T = the 
only non-zero term of the sum is the k = one (i.e., all the particles are condensed). 
From (16) one can immediately see that the critical value oiU/t for the superfluid- 
insulator transition is: Uc/t — Azn. 

2.2 Magnetism and Superfluidity 



We study now a simple model which displays coexistence and competition be- 
tween antiferromagnetism and Bose- Einstein condensation (superfluidity). The model 
represents a gas of interacting spin-1/2 Jordan- Wigner bosons with Hamiltonian 
{t>0) 

H = tY. (^Ua+H.c.)+j5^(si-Sj-^) + y^ninj-/2^nj, (17) 

(i,j>,o- (i,j> j 

where sj = ^\^^ iiv^jy is a s = | operator (<? denoting Pauli matrices). Notice that 
H is an extended t- J-like model of hard-core bosons instead of constrained fermions 
[10]. These hard-core bosons could represent three-state atoms, like the ones used in 
trapped Bose-Einstein condensates, moving in an optical lattice potential. For the 
sake of simplicity we will only consider the AF, J > 0, case. 

In the HL, H represents a Heisenberg-like Hamiltonian [9] in the presence of an 
external magnetic field {J^v = Jvn) 

H = Y,J^.S^^''(S)S''^{S) - l^'Y.^'\S) , (18) 

(i.j> j 

with Joo = y - J/2, Joi = Jo2 = t -hi = Ji2 = J22 = J/2, and ^' = §{2V - 
J/2) — /i. This HL furnishes the natural framework to analyze the symmetries of the 
Hamiltonian H. There is always an SU{2) spin symmetry generated by — 5^^, 
S^"^, and >S^^. When = and V = 2t, there are five additional generators of 
symmetries related to the charge degrees of freedom. Moreover, if J = V = 2t there 
is full SU{?>) symmetry. For 7^ 0, the only charge symmetry that remains is a 
U{1) symmetry generated by S^^ (conservation of the total charge). In this way 
the HL, leading to a unique OP from which all possible embedded orderings are 
derived, provides a unified description of the possible thermodynamic states of the 
system. Yet, it remains to establish the orderings that survive as a result of tuning 
the parameters of the Hamiltonian or external variables such as temperature and 
particle filling. 




Figure 1. Order fields A and i? as a function of tlie density po for different values of 
J/t and inverse temperature (3 = 10 (in units of The filled circle on the density 
axis indicates a quantum critical point. 



For arbitrary values of the parameters J/t, V/t, we do not know a priori how to 
determine exactly the phase diagram of H [12]. The resulting Hamiltonian {V = 2t 
with no loss of generality) is [10] 

^=-E(^4^ij +^ E -/^E^jo , (19) 

(i,j> cr^hi j 

where fx = zt — p, and the ordering operators 



bL = bLbl - bLbl 



(20) 

Vij "ia^jO "iO"jCT 

which transform as singlets with respect to the generators of SU{2) spin and charge 
symmetries, respectively: [tIJ., 5^2(21) (j)^^12(21)q)] ^ q ^ [^t^^^^^^ ^10(20) (i)^^io(20) 

The resulting MF Hamiltonian reads [10] 

H = -E[^(^u + ^ij) + E (^iij + ^-ij)] - E + ^ E 

(i,j> o-=T,i j i,a 

=E ^IcT^-k+Qi + ^^E ^L^-k+Qo + H.c. +(A - /x)nko +XY1 "ka],(21) 

where the sum of momenta k is performed over the reduced Brillouin zone (RBZ) 
with AF ordering wave vector Q, with = — 2JA7k, A^ = —AtB^-^, with 7k = 

\Yl,s'^^^^ nearest- neighbor vectors). Note that when 5 = in Hmf, the 

SU{2) spin and U{1), S^^, symmetries are conserved; the opposite case A = pre- 
serves 5^0*^0^^ + 5^o(°^) and S^^ + — S^^ symmetries. In Eq. (21) we have only 
considered homogeneous solutions [14]. 
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Figure 2. Schematics of tlie order parameter cliange (as a function of po) with the 
three su{2) axis describing 3 different su{2) subalgebras of 5^(3). Note that the plane 
t — I represents the 5 charge symmetry generators while the j | axis is associated 
to the remaining three generators of magnetism. 

The corresponding self-consistent MF equations to solve are 

J2 7k(&it^>lk+Qi)MF , 
^ keRBZ 

< Yl 7k(&L^-k+Qo)MF , (22) 

keRBZ 

V ^ keRBZ a 

We are thus left with a non-interacting system of SW bosons. Now we follow Colpa 
[15] and diagonalize para-unitarily the Hamiltonian matrix H. The application of a 
homogeneous linear transformation leads to [10] 

5 

H = ^ujiy, alk«*k , (23) 

keRBZ i=0 

where the mode energies are, at least, two-fold degenerate. In Fig. 1, we display 
the orders A and S as a function of po = 7^ X]j(^jo) at very low temperature and 
different J/t ratios for a two-dimensional lattice [16]. The relation between the OPs of 
the original problem, Eq. (17), and A and B is given by J2i j ^^^^^'~^^HHa^Sa') mf 
B"^, j^^^^^'~'^^H^t ^i) mf oc A^, and justifies the labeling of the phases AF 

(antiferromagnet) and BEC (Bose-Einstein condensate) in Fig. 1. 

A way to qualitatively understand this quantum phase diagram is to look into 
the OP space as displayed in Fig. 2. As po (or chemical potential) varies from to 1, 
the OP (depicted as an arrow) moves in OP space. When po = the order is purely 
AF and the arrow lies on the || axis. For po 7^ 0, the arrow has projections onto 
the 3 su{2) axis, i.e., the AF state coexists with a BEC state. There is a particular 
critical value of po = Poc < 1 for which the AF ordering vanishes and the OP is 



purely BEC with the arrow lying in the 10—^0 plane. This BEC ordering persists 
until po = 1) where it vanishes. 

There are some open issues. One regards the application of the HMFT approach 
to study fermionic problems, for example, a Hamiltonian like Eq. (17) but where the 
operators for different modes on a lattice i anticommute. The method could 
certainly be used, however, fermions do introduce a non-local gauge potential [9] 
leading to an effective dynamical frustration which is difficult to handle in a controlled 
manner. Another issue concerns the application of the HMFT method when longer- 
range interactions are involved. There is already evidence from work on the J1-J2 
SU{2) Heisenberg model [17] that our HMFT will work in those cases as far as 
homogeneous phases are concerned. Actually, the SW MF theory introduced by 
Arovas and Auerbach [18] in the fundamental representation is a particular case of 
our general HMFT. Finally, our methodology exhausts all broken symmetry instances 
but it is still quite possible to have purely topological quantum orders and their 
corresponding phase transitions which cannot be described by broken symmetries 
and associated OPs [19] and, thus, are not included in our framework. 

3. SUMMARY 

Summarizing, we developed a theoretical framework and a calculational scheme 
to study coexistence and competition of thermodynamic phases in strongly corre- 
lated matter. In our method (given a Hamiltonian modeling the physical system) 
the order parameters are not guessed but rigorously determined from group theoret- 
ical considerations as symmetry generators of a hierarchical language. In this way, 
the Hamiltonian operator (which does not necessarily have the full symmetry of the 
hierarchical group) is expressed in terms of symmetry generators. Then, in a non- 
phenomenological approach dubbed hierarchical mean-field theory, we approximated 
the dynamics (and thermodynamics) treating all possible local order parameters on 
an equal footing, i.e., without preferred symmetry axis. One could say that this 
procedure follows the guiding principles of maximum symmetry and minimum infor- 
mation. This allowed us to obtain in a simple manner the phase diagram of a model 
problem exhibiting a superfiuid-insulator quantum phase transition and another dis- 
playing coexistence and competition between antiferromagnetism and superfluidity. 
Combined with an analysis of the fluctuations (to analyze the stability of the mean- 
field) one now has a simple machinery to design phase diagrams. 

As can be inferred from our presentation, there are two complementary aspects 
to studying competition and coexistence between phase orderings in strongly coupled 
quantum systems. One is the direct discovery of the hidden unity and subsequent de- 
termination of the possible phases and their transitions, given a Hamiltonian operator 
modeling the complex material of interest. This is the aspect we have described in the 
present paper. The second aspect, to be discussed in a separate publication, involves 
the design or engineering of new states of matter (i.e., new quantum orders) using 
the inverse path of logic. Essentially, the idea consists of tailoring effective Hamil- 
tonians based upon a general symmetry analysis of the possible orderings one would 
like to realize at zero temperature. Tuning the parameters of these symmetry-based 
effective Hamiltonians allows one to move in parameter space along the previously es- 
tablished orderings. Indeed, this strategy finds its experimental realization in recent 



work done on atomic BEC systems in optical lattices [20] , and our approach provides 
a unique theoretical guidance to achieve that goal given the possibility of control and 
tunability of the interactions of the elementary constituents (i.e, quantum control). 

This work was sponsored by the US DOE. GO thanks the US Army Research 
Office for travel support. 
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